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Abstract. The moduli spaces of hyperbolic surfaces of genus g with n geo- 
desic boundary components are naturally symplectic manifolds. Mirzakhani 
proved that their volumes are polynomials in the lengths of the boundaries 
by computing the volumes recursively. In this paper we give new recursion 
relations between the volume polynomials. 



For L = (L±, L2, L n ), a sequence of non-negative numbers, let M g _ n (li) be 
the moduli space of connected oriented genus g hyperbolic surfaces with n labeled 
boundary components of lengths L±, L n . A cusp at a point naturally corresponds 
to a zero length boundary component. When L = 0, that is there are n cusps, the 
moduli space Ai gin (0) is naturally identified with the moduli space of conformal 
structures on a genus g oriented surface with n labeled points, also known as the 
moduli space of curves with n labeled points. The identification uses the fact that 
in any conformal class of metrics there is a unique complete hyperbolic metric, and 
for every conformal automorphism there is a corresponding isometry. 

On the moduli space .M 9jra (L) lives a natural symplectic form oj, defined precisely 
in Section [5] The volume of the moduli space is 



When (g,n) = (1, 1) we instead take half of the integral of u, an orbifold volume, 



which fits well with recursion relations between volumes, and relations with inter- 
section numbers on the moduli space. Mirzakhani uses the true volume of Mi i(Lx) 
in [HI El and includes an extra factor of a half in her formulae. 

Theorem 1 (Mirzakhani ©). V g . n (L) is a polynomial in L = (Li, L n ). The 
coefficient of L a = L ai ..L a " lies in 7r 6 s- 6 + 2r H a lQ, | a | = ai + ... + a n . 

Mirzakhani proved this using a recursion relation between volumes of moduli 
spaces: 



1. Introduction 





Mi.i(ii) 




(1) 



— (LxVg^L)) = A g , n (L) + B g , n (L) 
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where A g , n (L) consists of integral transforms of Vg_ x n+1 and B g , n (L) consists of 
integral transforms of V g . n ~i- The ' in V g '_ 1 n+1 allows disconnected surfaces which 
is an efficient way to encode the volumes of pairs of moduli spaces, or equivalcntly 
volumes of products M. gi: m x M g2: n 2 for gi + g 2 = g and n x + n 2 = n + 2. We 
have omitted the L dependence in V g l _ 1 n+1 and V^„_i because it requires further 
explanation. See Section l2~2l for precise definitions of A g . n (L) and 6 Sj „(L). 

The main idea of this paper is to use intermediary moduli spaces to give new 
recursion relations between volumes of moduli spaces. The intermediary moduli 
spaces consist of hyperbolic surfaces with a cone point of a specified angle. Hyper- 
bolic geometry is an ideal setting for studying cone points, although a cone point 
does make sense more generally in terms of a conformal structure on a Riemann 
surface. A cone angle of corresponds to a cusp marked point and as the cone 
angle goes from to 2ir this corresponds, in a sense, to removing the marked point. 
This leads to interesting relations between the moduli spaces. These intermediary 
moduli spaces are reminiscent of the moduli spaces of anti self dual connections 
with cone singularities around an embedded surface in a four-manifold, used by 
Kronheimer and Mrowka 0] to get relationships between intersection numbers on 
instanton moduli spaces. 

In it is shown that one can interpret a point with cone angle in terms of 
an imaginary length boundary component. Explicitly, a cone angle <fi appears by 
substituting the length i<f> in the volume polynomial. Mirzakhani's results, Theo- 
rem Hand (D use a generalised McShane formula (Hj on hyperbolic surfaces, which 
was adapted in JOj to allow a cone angle 4> that ends up appearing as a length 
icf> in such a formula, and hence in the volume polynomial. We do not describe 
the generalised McShane formula in this paper, although in Section l2~2l we give the 
underlying idea in terms of coordinates on the hyperbolic surface. 

Theorem 2. For L = (L%, L n ) 

(2) V g , n+1 {L;2m) =J2 L k V g , n (L)dL k 

k=i Jo 

and 

(3) ^f^(L,2Tri) = 2Tri(2g-2 + n)V g>n (L). 

We think of the theorem as describing the limit of the volume and its derivative 
when a cone angle tends to 2ir, and hence is removable, although the statement of 
the theorem is independent of this interpretation. The recursion relations J5J and 
(PI give information about the volume of M. g<n +x(Li, .... L n+ i) from the volume of 
the lower dimensional moduli space M. g ^ n (L\^ ...,L n ). This contrasts with Mirza- 
khani's relation (JTJ which uses many lower dimensional moduli spaces as described 
above. We discuss this difference more in Section0]and in particular give a simpler 
algorithm to determine Vo, n (L). 

There are three potential approaches to the proof of Theorem|21 one of which we 
present here, in Section |3 and the others we cannot yet do. Firstly, the theorem 
must follow from Mirzakhani's recursion relation Q since that relation uniquely 
determines the polynomials. Secondly, the theorem should follow from an analysis 
of the cone angle moduli spaces. Thirdly, the statement of the theorem is equivalent 
to relations between the coefficients of the volume polynomials which are intersec- 
tion numbers of ip classes and k classes (see Section for definitions), so relations 
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between the latter can be used to deduce the theorem. It is this third approach 
that we present here. 



2.1. Fenchel-Nielsen coordinates. Given an oriented hyperbolic surface £ with 
n geodesic boundary components, cut along a maximal set of disjoint embedded 
simple closed geodesies. The resulting pieces are genus zero hyperbolic surfaces each 
with 3 geodesic boundary components, i.e. hyperbolic pairs of pants. Each pair of 
pants contributes Euler characteristic -1, so there are 2g — 2 + n = — x(S) pairs of 
pants in the decomposition, and hence 3g — 3 + n closed geodesies (not counting the 
boundary geodesies.) We call such a decomposition a pair of pants decomposition 
of the surface E. A pair of pants decomposition gives a coordinate system on the 
moduli space of hyperbolic surfaces - the lengths of the disjoint embedded simple 
closed geodesies and angles which we will now define. 

Each boundary component of an oriented hyperbolic pair of pants has a pair of 
distinguished points - the closest points to the other two boundary components - 
which divide the boundary circle into two equal length pieces. A closed geodesic 
in a pair of pants decomposition contains two pairs of distinguished points coming 
from each side of the geodesic. The (oriented) angle between the pairs is a number 
9 e [0, 1], where I is the length of the geodesic, well-defined up to 9 i— > 9 + 1/2. 

So we have (h,9i) for i = 1, ...,3g — 3 + n, the Fenchel-Nielsen coordinates. 
Although 9i e [0, k) it is natural instead to lift it to 9i G K. The coordinates take 
values in 



which is known as Teichmiiller space. This statement also uses the fact that any 
triple of non-negative numbers give boundary lengths of a unique hyperbolic pair 
of pants and two oriented hyperbolic pairs of pants can be glued along boundaries 
of common length with any angle between distinguished points. 
The coordinates are canonical coordinates for a symplectic form 



which is a rather trivial statement. What is quite deep is the fact that the sym- 
plectic form is invariant under the action of the mapping class group Mod gn , of 
isotopy classes of orientation preserving diffeomorphisms of the surface that pre- 
serve boundary components, and hence descends to a quotient of Teichmiiller space 



known as the moduli space. The action can be understood by noticing that a 
given hyperbolic surface has infinitely many pair of pants decompositions that are 
geometrically different (but topologically the same and thus related by an element 
of Modg }Tl .) Each different decomposition gives different lengths and angles to the 
same hyperbolic surface, and hence different coordinates, whereas the symplectic 
form J3J is preserved. This gives rise to the volume of the moduli space Vg <n (L) 
given by the integral of uj 39 ~ 3+n / (3<? — 3 +ri)\ over the moduli space, or equivalently 
over a fundamental domain for Mod g ^ n in 7^„(L). 



2. Volume of the moduli space 



T S ,„(L) = (R+ x Rfa- 



(4) 




M 9in (L)=T gin (L)/Mod, 
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2.2. Coordinates on a hyperbolic surface. It is useful to view a hyperbolic 
surface from one boundary component chosen from the n boundary components. To 
motivate this idea, let us consider isometries of the surface that leave the boundary 
components invariant, although they may be rotated. A non-trivial isometry detects 
when the mapping class group does not act freely on 7^„(L). An important fact 
is that an isometry that fixes any boundary component is necessarily the identity. 
To see this, note that any point on the surface lies on a geodesic that meets the 
ith boundary component, say, at right angles - simply take the shortest path to 
that boundary component. (The path lies on the interior of the surface since if it 
was to meet another boundary component, it would meet at a positive angle - two 
geodesies cannot be tangent - and thus there would be a shorter path by taking a 
shortcut to the boundary near the angle.) An isometry that fixes the ith boundary 
component will preserve the orthogonal geodesic, so in particular it fixes all points 
along that geodesic. Applying this argument to any point on the surface shows that 
all points of the surface are fixed by the isometry. Thus any specified boundary 
component, say d\, detects non-trivial isometries. 

Another way to view the discussion in the previous paragraph is that the geodesic 
boundary component di gives a coordinate system on the surface - to any point 
on the surface assign its distance from d\ and the point on d\ where the shortest 
geodesic meets. More generally, take any geodesic beginning at a given point on 
the surface and meeting d\ perpendicularly, and assign to the point its length and 
the point it meets d\ . This makes the coordinate system locally smooth, at the 
cost of losing uniqueness for the coordinates of a point. 

Mirzakhani uses this coordinate system in the following way. Project points 
onto the second coordinate, which takes its values in d\. Now suppose that there is 
another boundary component, di say. The projection of di is an interval If C d\. 
More precisely, the projection is a collection of infinitely many disjoint intervals 
{I°i I 3 — 0; °°} since we take any perpendicular geodesic, not just the shortest 
one, resulting in non-unique coordinates. 

The sum of the lengths fi = • 1(1 f) is a well-defined function on the moduli 
space .A4 s>n (L). The length of a single interval l(I[) is well-defined on Teichmiiller 
space 7^„(L), and although it does not descend to the moduli space, l(Il) descends 
to an intermediate moduli space: 

^,n(L) 
_ I 

M g , n (L) 

I 

AWL) 

and Mirzakhani shows that this enables one to integrate the function fi = J^j 1(1 j) 
over J\4 g , n (L) yielding a polynomial, calculable from V g , n -\. The n—1 collections of 
intervals {I? \ j = 0, oo}, i = 2, n are disjoint from each other and Mirzakhani 
similarly shows that the complementary region (up to a measure zero set) gives a 
well-defined function f c on the moduli space which can be integrated in terms of 
lower volumes. Since f c + fi = the sum of all of the integrals gives 

/ Lidvol = L{V g , n (L) 

•Am 9 ,„(L) 
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the derivative of which can be calculated and leads to Mirzakhani's recursion rela- 
tion: 



— (LiV g ,„(L)) = A g , n (L) + B 9 ,„(L) 



d_ 
dLi 

For completeness we will define the right hand side although this will not be used 
further in the paper. Put L = {L2, ...,L n ) and let (L2, Lj, L n ) mean we 



remove Lj. Then 



Ag !n (L) = J K Ll (x,y)Vg_ 1>n+1 (x,y,L)dxdy 
where 

V g-hn+l( X >V>' L ) = V g-l,n+l(x,y,t) + V B1 ,™l ( X > L l) X V 92,n 2 (y, L 2 ) 

9i ,rii, Li 

and the sum is over all <?i + 52 = .9, "1 + "2 = n + 1 and Li U L2 = L. And 

n p 

#fl,n( L ) = ^ / K L l ,Lj(x)Vg 1 n-i(x,L 2 ,...,L j ,...,L n )dx. 
3=2 J 

The kernels are defined by 

K Ll (x,y) =H(x + y,Li), K Ll>Lj (x) = H(x, L x + Lj) + H (x, L x - Lj) 

for 

H{x,y)= l - f ' ' 



z V 1 + e 2 1 + e 2 , 
The derivation of these kernels comes from a detailed study of a hyperbolic pair of 
pants - the simplest hyperbolic surface to contain two boundary components. We 
refer the reader to 011] for full details. 

3. Characteristic classes of surface bundles 
3.1. Surface bundles. To any oriented topological surface bundle 

it i t Si i = l, ...,n 
B 

with n sections having disjoint images we can associate characteristic classes in 
H*(B), [5]. On X there is a complex line bundle 7 — > X with fibre at b G i? the 
vertical cotangent bundle T*7r _1 (6). A local trivialisation is obtained from a local 
trivialisation of the fibre bundle X. For each i — 1, ...,n pull back the line bundle 

7 to s*7 = 7,; — > -B. Define 

& = ci( 7i ) eH 2 (B). 

Let e = 01(7) G H 2 (X). (We use the terminology e because it is naturally the 
Euler class of 7. We have put a complex structure on 7 for convenience.) Define 
the Mumford-Morita-Miller classes 

k m =7ne m+1 G ff 2m (B) 

where m : H k (X) — > H k ~ 2 (B) is the umkehr map, or Gysin homomorphism, 
obtained by integrating along the (oriented) fibres. Alternatively, the umkehr map 
is obtained from the composition 

7H : H k (X) P S H d _ k (X) -% H d ^ k (B) -S H k - 2 {B) 
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where d — dimX and PD denotes Poincare duality. The Mumford-Morita-Miller 
classes ignore the n sections Si. Use instead the sequence 

tt, : H k c {X - U Si (B)) P S H d _ k (X - Ua^B)) ^ H d _ k (B) p -2 H k - 2 {B) 

where H k (X — Usj(S)) denotes cohomology with compact supports. Define the 
kappa classes 

K m =7r,e™ +1 G H 2m (B) 

where e c = e(^) 6 H 2 (X — UjSj(-B)) is the Euler class with compact support. It 
has the property that on any fibre E 

(e C) S - UiSi(B)) = -x(S - Us 4 (B)) 

which generalises (e,E) = — x(E). It is convenient to work with the compact 
manifold X and in place of e c use its image in H 2 (X) 

H 2 (X-UMB)) - # 2 (*) 

n 

e c ^ e„ = e + 53PD[« i (B)]. 

z=l 

The expression for e n is deduced from its two properties 

(5) (e n ,S) = -x(E-U i s i ( J B)), e„ • PD^iB)] = j = l...n 

the first because it is defined by restriction, and the second because it lies in the 
kernel of the map H 2 (X) -> H 2 (\J lSl {B)). 

We will need relations between classes obtained by simply forgetting a section. 
Now 

e„+i = e„ + PD[s n+ i(S)] 
so from e„+i • PD[g n+1 (B)] = and c = a + b^ c m+1 = a m+1 + &EJ1 cPa m ~ j 

e™# = e™ +x + PD[s n+x {B)]-a™ 
thus the forgetful map 7r„ + i induces 7r* +1 : H*(B) — > H*(B) satisfying 

(6) « m = n n+1 K m + ip™+i 
and the straightforward relation 

(7) =<+i^, J = 

To any E bundle n : X ^ B with n sections s$ corresponds the pull-back E 
bundle 7r*X — > X with n sections n*Si and a further tautological section s„+i. In 
some sense the section s n +i gives all possible ways to add an (n+l)st section to the 
bundle over B. In this context the forgetful map has two interpretations. As the 
map 7r* +1 : H*(B) — > H*(B) discussed above, and also as tt* : H*(B) — > H*(X). 
The two are related by 

S n+1 07T = K n +l- 

The pull-back relation JfjJ) looks the same for tt* 

(6a) K m = 7T*/t m -I- ip™ +1 

whereas the relation Q needs to be adjusted to 

(7a) ^ =n*4 !j +PD[ Sj (B)}, j = l,...,n. 
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We have yet to mention that the tautological section s n +i ■ A — > ir*X does 
not have disjoint image from the other sections Tt*Si. After blowing up to separate 
the images of the sections, we are naturally led to consider surface bundles n : 
X — > B that allow fibres with mild singularities. More precisely, the singular fibres 
may be stable curves - they consist of a collection of smooth components meeting 
at nodal singularities with the property that each component has multiplicity 1, 
and negative Euler characteristic after we subtract all labeled points and common 
points with other components. We call X a bundle of stable curves or simply a 
bundle with singular fibres, although strictly it is no longer a fibre bundle. The 
cohomology classes ipi an d n m extend to this situation. Their definitions are best 
understood when we put a continuous family of conformal structures on the fibres, 
or we assume the stronger property that X and B are complex analytic varieties. 
Define 7 = Kx®^* K^ 1 , essentially the vertical canonical bundle (relative dualising 
sheaf.) This coincides with the definition above on smooth fibres and generalises 
the definition to singular fibres. One can make sense of sections of this bundle along 
singular fibres in terms of meromorphic f-forms with simple poles and conditions 
on residues but we will not explain this here. The definitions of ipi and are 
as above. Relations Hfiajl and H7a(l generalise to bundles of stable curves. Proofs can 
be found in pQ and 

A simple example will demonstrate the definitions and relations. Let X be the 
blow-up of P 1 x P 1 at the three points (0, 0), (f , 1) and (00, 00). The map from X 
to the first P 1 factor realises AT as a surface bundle 



which we equip with four sections si(z) = (z, 0), S2(z) = (z, 1), Ss(z) = [z, 00) and 
Si(z) = (z, z). The general fibre is genus with 4 labeled points, and the singular 
fibres, at 0, f and 00, are stable curves with two irreducible components each with 
two labeled points (and a common intersection point.) We can generate ^(A) by 
H, F, Ei, E2 and E3 where Ei are the exceptional divisors of the blow-up and 
H = P 1 x {w} and F = {z} x P 1 for any w and z different from 0, 1 and 00. We 
use these curves to represent their divisor class, homology class and their Poincare 
dual cohomology class. Then 



Since X is the blow-up of the pull-back of the P 1 bundle over a point with three 
sections, Ij6a|) and (jTafl are also evident. 

3.2. Intersection numbers. Let us use M. g , n to notate the moduli space of genus 
g curves with n labeled points, which is isomorphic to the moduli space of genus 
g hyperbolic surfaces with n labeled cusps, A4 g , n (L) with L = 0, and M g , n the 
Deligne-Mumford compactification which adds stable curves to M g , n - Wolpert [T3] 
showed that the symplectic structure to on M g , n extends to M.g t n- The ipi and 
n m classes naturally live in H*(M g , n )- They are associated to a universal surface 
bundle over Mg, n , essentially given by M. gi n+i with map forgetting the last labeled 



A 

pi 



■A 
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point, and any bundle X equipped with conformal structures on fibres is the pull- 
back of the universal bundle under a map B — ► M g , n - 

Theorem 3 (Mirzakhani). The coefficient C a of L 2 " 1 ...L^ n in V g>n (L) is 

( 8 ) c « = ^nro — it, itttt / ^-C^ 3s - 3+ ^ |a| 



2l«la!(3 5 - 3 



M„ 



This is proven in by showing that A4 g ,n(L) is the symplectic quotient of a 
larger symplectic manifold by a Hamiltonian T n action, where a fixed value of the 
moment map corresponds to fixing the lengths L±, L n of the geodesic boundary 
components. Any such quotient is equipped with n line bundles coming from the 
T n action, and their Chern classes are related to the coefficients of the volume 
polynomial. In [7] Mirzakhani used this together with her recursion relation for 
the volume polynomials to give a new proof of Witten's conjecture [IT] regarding 
intersections of ip classes on M. 9in . In the original proof of Witten's conjecture, 
Kontsevich [3j calculated the Laplace transform of the top degree terms of V^ n (L). 
It would be interesting to understand the Laplace transform of the whole polynomial 

In the following, write tp a for tp" 1 ...ip" 71 and ignore the term if there is an ctj < 0. 
For ease of reading, note that in all formulae the variable j sums from to to while 
the variable k sums from 1 to n. 



Lemma 1. The equation 

(L,27ri) = > / L k V g , n (L)dL 
t_i Jo 



k=l ■ 



is equivalent to 



(9) E(-!) J m /- «+i<" i = E/_ < i -C fc - 1 -c-< 

„•_„ \3 J JMrr.rr + l t_1 J Man 



J=0 ^ J ' "J^g.n+l k=l • 

for all a and to. 

Proof. Assume that \a\ + to = 3g — 2 + n since otherwise (0 is zero on both 
in Vg tn+ i(Ij,2iri) is 



sides. By © and substitution of L I 2 / +1 with (27ri) 2j , the coefficient of L 2ai r2a " 



2I q I+Jq;!7'!(to - 



3 



2l Q l+Ja!j!(TO - j)! Jm s ,„ +1 



nra — |a|_2m m 

2 ^ Vr-ivf'") / ,h«,ij..*™-i 



y\ to! ^— ' 



where we have used the identity cj = 2it 2 ki proven in |12| . 
The coefficient of L 2ai ...L 2 n a ~ in L k V g , n dL k is 

^ / I Oi\ 1 /Ctrl 7i 



2l«l- 1 (2a fe )a!TO! ^ 
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a\ to 



Add this expression over k = 1, n and divide both sides by the factor 
2 m ~\ a \ir 2m / a\m\ to prove the lemma. □ 



Lemma 2. The equation 
BV 

a 5'" +1 (L, 2m) = 2m(2g - 2 + n)V g , n (L) 



< + l 

is equivalent to 



(10) L rtf&*r j = (2 9 - 2 + n) / 

? =0 \ J / JMg, n +l J A 



(2 5 -2 + n) /_ ip a K? 

M„. n 



Proof. The proof is much like the proof of the previous lemma. The coefficient of 
Ll a \..L 2 n a « in dV g , n+1 /dL n+1 (L,27ri) is 

2 l-l+i+ia!(j + l)!(m - j)! i^,„ +1 ^ 

om— |a| 2m m / m \ /• 

and the coefficient of L\ ai ...L^ n in V g , n is 



'9,' 

a! m! 



M g , n 



so the equivalence follows. 



□ 



Completion of the proof of Theorem It suffices to prove the relations (|5J) and 
(|f Op t . Notice that when to = 0, © and (|I0|l are respectively the string and dilaton 
equations which were proven by Witten in JT] . The method of proof for the more 
general identities is similar. 

Let 7r : X — > B be a bundle of stable curves with n disjoint sections and tt*X 
the pull-back bundle with n + 1 sections. Blow up 7r*JT along the intersections 
of images of sections to get a bundle of stable curves over X with n + 1 disjoint 
sections si, s n _|_i. Our aim is to compare ?/> and k classes in H*(X) and H*(B). 

Take the integrand of the left hand side of © and consider its image under the 
umkehr map. 

"i { p-tf (7) n <c } = "t { c«i - n c fc } 



7Tl 



k=l 



To get from the first line to the second line we have used the pull-back formulae H6a|) 
and (7a|) and the fact that 7r* is a ring homomorphism, so in particular (71*77)™ — 
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7r*(?7 m ). To get from the second line to the third line we have used the fact that 
7T! : H*(X) — > H*(B) is an H*(B) module homomorphism, i.e. m{£,^*r{) = n\(£)r], 
together with the explicit evaluations 

7n(l)=0, 7r ! (s i (B)) = l 

most easily calculated from the Poincare duality description of m. Thus, in the 
product the image under 7n of the highest degree term 7t*(k™V' q ) is zero, the image 
of the second highest degree term constitutes the expression in the third line, and 
the lower order terms vanish since they contain products PD[sj (B)]-PD[sk(B)] = 
because the images of Sj and Sk are disjoint. 
Since 



T) = / 7n?7 

X JB 



choose X = Mg, n +i and B = M g , n , so © follows. 

The proof of t| 1 ( )f) is similar. Again apply the umkehr map to the integrand of 
the left hand side of (p3J). 

m / \ n j ( n \ 

D-i) j rw + \<- 3 n ^ k = - • (-1 - ^ + ir n r k * 

3=0 ^ ■' ' k=l J I k=l J 

= tt! |^ l+ i • (tt^D f[ (**1>k h + PD[s k (B)] ■ tt*^ 1 )} 

n 

= (2.g-2 + n)< i n^ fc - 
fc=i 

To go from the second line to the third line note that VVi+i coincides with the 
twisted Euler class e„ + i that satisfies (0 and hence 7ri7/; n+1 = 2g — 2 + n and 
^n+i • -P-D[sfc(-B)] = 0. Thus the top degree term constitutes the expression in the 
third line, and all lower degree terms vanish. □ 

Equations (J2J) and suggest that a direct analysis of the moduli space of cone 
surfaces with cone angle 9 w 27r, or more accurately an infinitesimal analysis near 
27T, will gives rise to intriguing phenomena. Equation ||3J) seems plausible since 
the removed cone point is free to wander around each hyperbolic surface with area 
27r(2g — 2 + n), so the change in volume is related to integrating over the smaller 
moduli space and along each fibre. Intuition for equation (|5J) seems less obvious. 



4. Use of recursion relations 

4.1. Classical volumes and low genus algorithms. The relations J5J and (J3J) 

can be used to give convenient information about the volumes of moduli spaces. 
In the following, when there is only one boundary component we use the variable 
L = Li. 

(i) When there is exactly one marked point the volume factorises 
(11) V g ,i(L) = (L 2 + 4tt 2 )P s (L) 
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(ii) The classical volumes of moduli spaces, i.e. when there are no marked 
points, are encoded in Mirzakhani's volume polynomials. In terms of the 
polynomial P g defined by 

(12) V„ -*•<*>- 1 



L=2-iri 



9 - 1 (.9 - 1) L 2 + 4tt 2 
(hi) The relation uniquely determines Vo, n +i from Vo n - 
(iv) The relations @ and © uniquely determine Vi n+i from V\ tn . 
Equations and (|12[) follow from the proof of Theorem [21 In that proof we see 
that when n = the right hand side of J5J) must be zero, which can be interpreted as 
there being no boundary lengths to integrate, and hence V 9t i(L) possesses a factor 
of L 2 + An 2 . To get (JT2J), use © together with V g ,i(2m) = 0: 

V 9 AL) 



2m(2g - 2)V gfi = ^ 



L=2iri 



lim 

L^i-Ki L — 2m 



L^i-Ki L 2 + An 2 a 

Proof of (Hi). This follows from the elementary fact that a symmetric polynomial 
f(x±, x n ) of degree less than n is uniquely determined by evaluation of one 
variable at any a S C, f{x±, ...,x n -i, o). To see this, suppose otherwise. Any 
symmetric x„) of degree less than n that evaluates at a as / does, satisfies 

f(xi,...,x n -i,a) - g(xi,...,x n -i,a) = (x n - a)P{xx, x n ) 

Th 

= Q(xi,...,x„) Y[(xj -a) 

3 = 1 

but the degree is less than n so the difference is identically 0. 

The volume Vo,n+i is a symmetric degree n — 2 polynomial in L 2 , .... so it 
is uniquely determined by evaluation at L n+ i — 2m\ and this is determined by Vo,n 
via ©. □ 

The proof of (iv) is similar to the proof of (hi). The degree of Vi >n +i is equal to 
n + 1 so the proof of (iii) shows that (2) determines Vi_ n+ i from Vi_ n up to the con- 
stant c in Vi, n+ i+cJ|" =1 (L^ + 47r 2 ). Now use (0 to determine c, and hence Vi,„ + i. 

Statement (iii) can be converted to an algorithm for calculating Vo, n - The al- 
gorithm using (J2J) turns out to be much more efficient than the algorithm coming 
from Mirzakhani's relation Q in genus 0, which needs Vb, n -i and pairs Vb,m; 1 l / o,n 2 
for all n\ +ri2 = n+ 1, to produce Vo,n- We have included a simple MAPLE routine 
in the appendix for calculating Vb, ra using J5J. (The notion of a "more efficient" 
algorithm is not so precise here. We have merely compared the speeds of different 
calculations on MAPLE.) 

In genus 0, the string equation - with m — - leads to an explicit formula for 
the top coefficients, or equivalently the following formula for genus intersection 
numbers without kappa classes: 

/ ^. = 3 

JM g , n \a>i,...,a. 

It seems reasonable to guess that when g = the relation © leads to an explicit 
combinatorial description of all genus intersection numbers with powers of ki , or 
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equivalently all coefficients of Vo,n ■ Zograf nas recursion relations between the 
constant coefficients Vb, n (0). 

4.2. Higher derivatives. We expect to have expressions for higher derivatives 
d k Vg in +i/dLj l+1 evaluated at L n+ \ = 2iri. Evidence comes from the fact that (J2J 
and (J3J) use generalised versions of the string and dilaton equations. The Virasoro 
relations are a sequence of relations for the top degree terms of V^ >n (L), with 
first two relations in the sequence the string and dilaton equations, so may also 
have versions in terms of evaluations of derivatives of the volume polynomial at 
L n+ \ = 2tti. The Virasoro relations recursively determine the top degree coefficients 
of the volume polynomials by using the relations in a clever way. In recent work 
UJ, Mulase and Safnuk showed how to extend the Virasoro relations to the full 
volume polynomials. It would be desirable to instead determine the polynomials 
recursively by relying on the more familiar fact that the derivatives of a function 
evaluated at a point determine the function. It would be interesting to know if 
one can express the results [5] in terms of derivatives of the volume polynomial at 
L n+ i = 2iri. 

In principle, we can use Mirzakhani's recursion relation to get expressions for 
higher derivatives of the volume evaluated at L n+ i = 2-ki. Differentiate the equation 

d(L n+ iVg^ n +i) 



dL r 

to get 



°Ag, n+1 ~l~ Bg, n+1- 



9 (£n+lK),«+l) _ dAg, n+1 <9g 3 ,n+l 

9L^ + i dL n+ i dL n+ i 

and evaluate at L n+ \. Substitute the equation for the first derivative, to get the 
following equation for the second derivative. Put £ = Yl^i^jd/dLj, the Euler 
vector field: 



d 2 v c 



^±i(L, 2m) = £ ■ V g . n (L) - 2ni(Ag - 4 + n)V g , n (L). 



By taking higher derivatives of Mirzakhani's relation we can recursively get equa- 
tions for higher derivatives. The strength of and © is the simplification of 
Mirzakhani's relations Q). It is not clear that the higher derivative relations ob- 
tained by the method above possess this same strength. 

Appendix 

maple routine for calculating vb, n - 

> # input: symmetric polynomial f in n variables Ll,...,Ln 

# output: symmetric polynomial S in n+1 variables LI , . . . ,L(n+l) 

# satisfying S(L(n+l)=0)=f 

sym:=proc(f) local i, j ,k,m,S, T,T1, prod, sum, epsilon: 
S:=f: 

epsilon : =array [1 . . 100] : 

for i from 1 to 100 do epsilon [i] : =0 : od: 
while epsilon [n+1] <1 do 

T:=subs(seq(L| I j=(l-epsilon[j] )*L| I j , j=l . .n) ,f ) : 
T1:=0: 

for i from 1 to n do 
prod:=l : 

for j from i+1 to n+1 do 
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prod:=prod*(l-epsilon[j] ) 
od: 

Tl:=Tl+prod*subs(L| I i=L I I (n+1) ,T) : 
od: 

sum:=0: for k from 1 to n do sum:=sum+epsilon[k] od: 

S:=S+(-l)~sum*Tl: 

for k from 1 to 100 do 

if epsilon[k]=l then epsilon [k] : =0 

else epsilon [k] : =1 : k:=100 end if: 

od: 

od: 

S:=simplify(S) : 
end: 

> # calculate the genus zero volumes recursively from evaluation 
# of V_(0,n+1) at L(n+l)=2*Pi*I 
for n from 3 to 12 do 
P:=0: 

for j from 1 to n do 
P:=P+int(L| I j *V [n] , L | | j ) 
od: 

Q0:=P: 

CO:=simplify(coeff (Q0,Pi,0)) : 

sim: =sym(C0) : 

V[n+1] :=sim: 

for k from 1 to n-2 do 

P| |k:=sim-C| I (k-1) : 

Q| |k:=subs(L| I (n+l)=2*Pi*I ,Q | | (k-l)-Pl I k*Pi~ (2*k-2) ) : 

C| |k:=simplify(coeff (Q| |k,Pi,2*k)) : 

sim: =sym(C Ilk): 

V[n+1] :=V[n+l]+sim*Pi~(2*k) : 

od: 

od: 
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